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ABSTRACT 



A factorization algorithm for a patron shower model based on the evolution of 
momentum distributions proposed in a previous work is studied. The scaling vio- 
lation of initial state parton distributions is generated using parton showers to an 
accuracy of the next-to-leading logarithmic (NLL) order of quantum chromodynam- 
ics (QCD) using the information from only splitting functions and initial parton 
distributions at some fixed low energy. In the algorithm proposed in this paper, 
the total momentum of the initial state partons is conserved in any factorization 
scheme for mass singularities. As an example, the scaling violation of the parton 
distributions and the transverse momentum distributions due to initial state parton 
radiation are calculated with the jet calculs scheme. 

§1. Introduction 

In perturbative quantum chromo dynamics (QCD), logarithmic contributions due 
to collinear parton production are subtracted from a hard scattering cross section, 
and these contributions are absorbed into the parton distributions of a hadron. In 
higher order calculations, the factorized terms from the hard scattering cross section 
are compensated for by associated terms of the splitting functions. Therefore, the 
scaling violation of the parton distribution functions also depends on the factorization 
scheme. The change of the hard scattering cross section is compensated for by the 
change of the distribution functions. We study this situation using a Monte Carlo 
calculation. 

In a previous work, 1 ) an algorithm that can reproduce the scaling violation for 
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flavor singlet partons is constructed to the next-to-leading logarithmic (NLL) order 
of QCD. That algorithm consists of a model based on the evolution of momentum 
distributions. In this model, a scaling violation of the parton distributions is gener- 
ated by using only information from the splitting functions of the parton branching 
vertices and input distributions at a given energy, Qq. It has been found that, in 
the MS scheme, 2 ^ the method reproduces the scaling violation of the flavor singlet 
parton distributions, up to their normalizations. 

In the model, the total momentum of the initial state partons is conserved. 
Furthermore, it is not necessary to introduce non-trivial weight factors in order to 
reproduce the energy scale dependence of the momentum distribution of partons in 
the initial state. These are advantages of this type of model in the calculation of the 
parton evolution to the NLL order of QCD. 

The purpose of this paper is to extend the parton shower model proposed in 
Ref. 1) in order to generate parton evolution with any factorization scheme. In §2, 
the algorithm for the mass singularity factorization of the parton shower model is 
explained. As an example, the explicit expressions for the jet-calculus scheme are 
presented in §3. In §4, some numerical results obtained using the proposed method 
are presented. Section 5 contains a summary. 



The parton evolution depends on the subtraction scheme used for the mass 
singularity, which appears due to collinear parton production. In the D = 4 — 2e 
dimensional integration of the phase space, the factorized contribution in 0(a s ) with 
a factorization scheme F can be written 



where i and j indicate the types of partons (q for a flavor singlet quark and g for 
a gluon). Here, M denotes the factorization scale of the collinear contributions. 
In 4 — 2e dimensions, the strong coupling constant is defiend by a s fi 2e , with the 
dimensionless coupling a s . The mass scale [i is also identified with the renormal- 
ization scale of the strong coupling constant. Here, z is the momentum fraction of 
the daughter parton and 7^ is the Euler constant. Finally, P-j\z) represents the 
infra-red unregulated splitting functions at the leading logarithmic (LL) order of 

QCD. 3 ) For example, AF^ S] (z) = 0, with fi 2 = M 2 , in Eq. (2-1) corresponds to the 
factorized term with the MS subtraction scheme. 2 ) 

In order to preserve the total momentum of the initial state partons, we subtract 
the infra-red regulated terms 



§2. Collinear parton factorization 



FiP(e,z,M 2 /, 2 ) = F^(e,z) + /f (.)log (^) (.), (2-1) 



with 




(2-2) 



FW(e,z,M 2 / f , 2 )= 1 -\zFl?(e,z,M 2 /v 2 ) 



+ 
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6(1 -z) I' dyypW(e,z,M 2 /» 2 ), (2-3) 
J o 



F^(e,z,M 2 /^) = F^(e,z,M 2 /^ 2 ) 
F^(e,z,M 2 / f , 2 ) = F^(e,z,M 2 / f i 2 ) 



(2-4) 
(2-5) 



F^(e,z,M 2 /^ 



F^(e,z,M 2 /^ 2 ) 



- 6(1 - z) f 1 dyy2N f FVp (e, z, M 2 //x 2 ) (2-6) 

JO 

from the hard scattering cross section. Here Nf denotes the number of flavors, and 

(2-7) 

for an unregulated function f(z) at z = 1. ^.Frorn Eqs. (2.3)-(2.7), we find that 

/ 1 ^[Fjfl(e,z,M 2 /M 2 ) + i ; J?(€,«,M 2 /M 2 )l =0 (2-8) 
Jo L J 



f(z) = f(z) - 6(1 - z) j dyf(y) 

o 



and 



f 1 dzz \2N f F$ ( e , z , M 2 /fi 2 ) + Fjjl (e, z, M 2 / ^ 
Jo L 







(2-9) 



are satisfied for any factorization scheme F. 

Although M and fi are independent parameters, we set M = ji for simplicity in 
the following calculation. 

Including the NLL order terms, the splitting functions, P\, (a s , z), are given by 



a s (M 
2vr 



2\\ 2 



(2-10) 



where P^f \z) and p\J ](1) (z) are the infra-red regulated splitting functions of the LL 
order and of the NLL order, respectively, calculated with a factorization scheme F. 
Using 



AFlP(z)=FlP(e,z,l)-F^(e,z) 



with e — > 0, the splitting functions to the NLL order are given by 



P^\z) = Pf^\z)-%AFf(z). 



(2-11) 



(2-12) 



Here, fio = 11 — 2/3Nf, and P^ AS ^ 1 \z) is the splitting function calculated with the 

MS scheme. 4) 

The quantities 



pW(a s ,z)=pVp(a a ,z) + pW{a,,z) 
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(2-13) 



4 



and 

Pf\a s ,z) = 2N f pW(a s ,z) + pW(a s ,z) (2-14) 

satisfy 



dzzP\ F] {a s ,z) = (2-15) 



for i = q,g. This guarantees that the total momentum of the initial state hadron is 
conserved, that is, 

f dxx [U F (x,K 2 ) + G F (x,K 2 )] = 1 (2-16) 
Jo 

for any K 2 , where x is the momentum fraction of a parton inside the hadron. Here, 
Sf{x, K 2 ) and G F (x, K 2 ) are the particle number distribution function of the flavor 
singlet quarks and that of the gluons in the factorization scheme F, respectively. 
Using 

f 1 dzzP\ F \z) = - f 1 6 ' dzzP\ F] {z), (2-17) 
Jl-S Jo 

the non-branching probablities for the flavor singlet quarks (i = q) and for the gluons 
(i = g) are defined by 



n^(K 2 ,K 2 )=e W 



f Kl dK 2 .1-6 



(2-18) 



Here, the partons inside an initial state hadron have spacedike virtualities (i.e., 
kf = —K 2 < 0), and 5 denotes the resolution of the momentum fraction of the final 
state partons. Although in the Monte Carlo calculation, only unregulated splitting 
functions are used, the infra-red regulated terms presented in Eqs. (2.3)-(2.6) must 
be subtracted from the hard scattering cross section in order to cancel the dependence 
on the factorization scheme in the hadronic cross section. 

Using this algorithm, the momentum distributions are reproduced up to their 
normalizations without the introduction of any non-trivial weight factor. Therefore, 
the model guarantees the total normalization of the momentum distributions for the 
flavor singlet partons in any factorization scheme. 

§3. Factorization in the jet calculus scheme 

In the jet calculus (JC), 5 ) the singular terms due to the collinear parton radi- 
ation, which are integrated over the t-channel momentum of the radiated parton, 
with the range of integration < — t < M 2 , are given by 



F\-\e,z,M 2 /^) 



r(i - e) 



M 2 



— e 



(3-1) 
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with 



(3-2) 



for i,j = q,g. Here, the functions P[j(z) are given in Ref. 6). 



From Eq. (2-11), the functions Af\^ C \z) in the JC scheme are given by 



Kl 



AF\f\z) = log(l -z)- P' qq {z) 

AF\¥\z) = P^(z)\og{l -z)- P' gq (z), 
^ g C \z)=P£\z)\og{l-z)-P> qg 



AF q f\z)=P£\z)\og{l-z)-P' qa (zl 
AFlf\z) = - z [zP^(z)\og(l-z) 



+ -S(l-z)-C F , 
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+ S(l-z)2N f T R -, 



with C F = 4/3 and T R = 1/2. 

In the Monte Carlo calculation, we study the following two cases. 



(Case 1): The splitting functions 



2vr 



with 



p[JCi](i)^ = p[MS](i) 



2vr 



(*)-|af, 



^ci](i) Wj 



[JC] W, 



(3-3) 

(3-4) 
(3-5) 

(3-6) 



(3-7) 



(3-8) 



are used in the non-branching probabilities, and the momentum fraction z is gen- 
erated according to the probability zP^ C1 ^ (a s , z). 



(Case 2): AF qq C ^ and APfg^ behave as ~ log(l — z)/(l — z) for z — ► 1. This may 
yield large higher-order contributions for z ~ 1. Such dangerous terms can be 
absorbed into the strong coupling constant 5 ) as 

a s (K 2 ) - a s " ^)^ 2 ) , 



with 



fij(z) = exp 



The splitting functions in the JC scheme are modified as 



4/ C2] (« s ,,) 



2vr 



, <*8 (K'tj 
-P^(z)+ l 



2vr 



(3-9) 
(3-10) 



(3-11) 



with K^. = fij(z)(l — z)K 2 ~ (0.4 ~ where p F is the transverse momentum 

of the generated parton with space-like virtuality. The factors fij(z) are presented 
in Fig. 1. 
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Fig. 1. The factors /^(z) defined in Eq. (3.10). 



§4. Numerical results for the JC scheme 



In this section, some numerical results obtained using the algorithm explained 
in the previous sections are presented. 

In order to check the scaling violation generated by the Monte Carlo calculation 
with the JC scheme, we calculate the momentum distribution functions in an 0(a s ) 
approximation. The momentum distribution functions in the JC scheme are given 
by 



xU JC (x, M 2 ) = xE m (x, M 2 ) + 



a s (M 2 
2vr 



dz 



AFf^z^M 2 ) 



+ AF\f\ z) *G m (^M 2 ) 



(4-1) 



and 



a s (M 2 



2vr 



xG JC (x, M ) = xG m (x, M ) + 

+ AFlf\z)*G m (?- z ,M< 



dz 



AFlf\z)^ m ^M 2 ) 



(4-2) 



The distribution functions xEj^(x, M 2 ) and xGy^(x, M 2 ) are obtained with GRV(98). 7 ^ 
In the Monte Carlo calculation, the parton showers start from the momentum 

distributions at Mq = 5 GeV 2 . In Fig. 2, the momentum distribution functions for 

the gluon at M 2 = 10 3 GeV 2 are indicated by '□' for case 1 and by '+' for case 2 in 

the JC scheme, with 5 = lljK 2 (ll = 0.2 GeV 2 ). 

Here, the results obtained from the Monte Carlo calculation are plotted as 

dN/dyF, with yp = — lnx^, where dN denotes the average number of events gener- 



7 



ated within the range uf ± The quantity dN divided by dyF corresponds to 

the momentum distribution multiplied by xf [i.e., xf(xf/{xf, M 2 ))]. 

In the figure, the dotted curve and the dashed curve are the results obtained with 
GRV(98) in the MS scheme for M 2 = 5 GeV 2 and M 2 = 10 3 GeV 2 , respectively. The 
dash-dotted curve and the solid curve denote the distribution functions calculated 
using Eq. (4.2) for M 2 = 5 GeV 2 and M 2 = 10 3 GeV 2 , respectively. In Fig. 3, 
similar results for the flavor singlet quarks are also presented. 
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Fig. 2. The momentum distribution functions (multiplied by xf) for the gluon obtained with the 
algorithm explained in §§2 and 3. The results obtained from the Monte Carlo simulation are 
indicated by '□' for case 1 and '+' for case 2 in the JC scheme with Iq = 0.2 GeV 2 at M 2 — 
10 3 GeV 2 . The solid curve and dashed curve represent the corresponding results for the JC 
scheme and for the MS scheme obtained from GRV(98). The input distributions at Mq = 5 GeV 2 
are represented by the dotted curve for the MS scheme and the dash-dotted curve for the JC 
scheme. 

Next, we present the transverse momentum distributions of the generated par- 
tons in the JC scheme. The transverse momenta for the partons are restricted by 
the NLL order terms. ^ At the LL order of QCD, the upper limit of the virtuality 
for the generated parton is restricted only by a large energy scale, which contributes 
to the branching vertex. In the actual Monte Carlo generation, the allowed kine- 
matic boundary for the two-body branching is determined by the condition that the 
squared transverse momenta of the generated partons be positive. 

In Ref. 1), the transverse momenta of the initial state partons with space- 
like virtuality are generated using the effective two body vertices, where the three- 
body decay functions calculated in Refs. 8) and 9)*) are included as the boundary 
conditions of the virtuality for the out-going partons. 

The transverse momentum of the parton with momentum k% for the branching 
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*' The function F(x) in the appendix of Ref. 9) is given by F(x) = x 2 + (1 - x) 2 . 
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Fig. 3. The momentum distribution functions (multiplied by Xf) for the flavor singlet quarks with 
the algorithm explained in §§2 and 3. The notation is the same as in Fig. 2. 

a(p) — ► b(k\ + A^) + c{kz) is given by 



X3 is the momentum fraction of the parton with momentum A^. The virtuality K 2 A 
is restricted by the condition K\ < /a{—s) at the NLL order, instead of K\ < 
Vs/zsi—s), by the kinematic boundary of the two-body branching, due to the fact 
that p\ > for — p 2 <C — s,K\. Here, the factor of Ja for the two-gluon radiation 
processes is given in Ref. 1). 

In Fig. 4, the transverse momentum distributions for the gluons with final 
momentum fractions satisfying 0.09 < xf < 0.11 are presented for M 2 = 10 3 GeV 2 . 
The results obtained from the Monte Carlo simulation are represented by '+' for 
case 1, 'x' for case 2 in the JC scheme, and '□' for the MS scheme, with Iq < K 2 A < 
/a(— s). The solid curve and the dashed curve correspond to the case Iq < K 2 A < 
ys/x^i—s) due to p\ > and the parton radiations with fixed K\ = Iq (= 0.5 GeV 2 ) 
at the LL order, respectively. In Fig. 5, similar results for the flavor singlet quarks 
are presented. *) 



In this paper, we studied a factorization algorithm for the parton shower model 
based on the evolution of the momentum distributions, which was proposed in a 



The transverse momentum distributions with Iq = 1 GeV 2 are presented in Figs. 6 and 7 of 
Ref. 1). 



where K\ = 




§5. Summary 
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Fig. 4. The transverse momentum distributions for the gluons with final momentum fractions sat- 
isfying 0.09 < x F < 0.11 at M 2 = 10 3 GeV 2 . The results obtained from the Monte Carlo simu- 
lation are represented by '+' for case 1, 'x' for case 2 in the JC scheme with I 2 , < K\ < /a(—s), 
and '□' for the MS scheme. The solid curve and the dashed curve represent the phase space con- 
dition due to the relation p% > and the radiation of partons with fixed K 2 A = 1% (= 0.5 GeV 2 ) 
in the LL order of QCD, respectively. 
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Fig. 5. The transverse momentum distributions for the flavor singlet quarks with final momentum 
fractions satisfying 0.09 < xf < 0.11 at M 2 = 10 3 GeV 2 . The notation is the same as in Fig. 4. 



previous paper. In this algorithm, we subtracted the singular terms due to the 
collinear parton production from the hard scattering cross section so as to preserve 
the total momentum of the initial state partons with any choice of the factorization 
scheme. 

The results obtained with the jet calculus scheme were compared with those 
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obtained using the corresponding distribution functions calculated in an 0(a s ) ap- 
proximation. These two methods give consistent results for the scaling violation 
of the flavor singlet parton distributions. The algorithm developed in this paper 
reproduces the dependence on the factorization scheme of the parton distribution 
functions. 

For the transverse momentum distributions, in the large px region, we find that 
the dependence on the factorization scheme at the NLL order is smaller than the 
ambiguity resulting from the choice of the kinematic boundary at the LL order of 
QCD. 

The method presented in this paper may be useful in the construction of Monte 
Carlo generators beyond the LL order of QCD, which are desired for future high 
energy experiments, such as LHC at CERN. 
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